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Abstract. We study a three-component superfluid Fermi gas in a spherically 
symmetric harmonic trap using the Bogoliubov-deGennes method. We predict a 
coexistence phase in which two pairing field order parameters are simultaneously 
nonzero, in stark contrast to studies performed for trapped gases using local density 
approximation. We also discuss the role of atom number conservation in the context 
of a homogeneous system. 
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1. Introduction 

Multicomponent ultracold Fermi gases allow the study of several interesting questions in 
many-body quantum physics. In particular, understanding three-component pairing can 
reveal some properties of multi- or two-component pairing. In three-component systems 
the pairing energy does not only compete with temperature effects or Fermi surface 
mismatch energy but also with other pairing gaps. As in imbalanced two-component 
systems, non-BCS pairing mechanisms such as Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) 
[HE], breached pairing (BP) [3 J and phase separation phases are expected. 

Ultracold Fermi gases have opened up a way to explore multi-component gases 
experimentally. Recently, a degenerate three-component gas was successfully created [U 
[5]. The stability of three-component gases is hindered by the three-body recombination, 
reducing the pairing in the gas and the lifetime of the sample [Hid IB]- However, there are 
ways to stabilize the gas against such losses using for example optical lattices. Optical 
lattices are interesting also due to the rich phase diagram: theoretical investigations 
have found that color superconductivity competes with normal phase and formation of 

trions p eh eeu na na m ng. 

Both the SU(3) symmetric model, in which the different components have identical 
properties, and the non-SU(3) symmetric case have been widely studied [TH1 CEO HB1 
H91 EQl EH E21 E31 E3 ES]. The particularly important special case of SU(3) symmetry 
can be realized using alkaline earth atoms [26] or in optical lattices. However, with 
alkaline earth atoms the numbers of atoms in different components are not well defined, 
and only the total number of atoms is conserved. In contrast, the hyperfine energy 
spacing in alkaline atoms stabilizes the atom numbers, making the atom number of each 
component separately conserved. This is the case of most two-component Fermi gas 
experiments and a natural extension of these studies is a non-SU(3) symmetric case in 
which the atom numbers in all three components are fixed. 

Here we study such a three-component system with fixed atom numbers in all three 
components in a spherically symmetric harmonic trap using the Bogoliubov-deGennes 
(BdG) equations. We study the coexistence of the pairing gaps in these systems and 
discuss the scaling of the system size up to the thermodynamical limit. 

In section [2] we give an overview of the three-component system and the 
corresponding mean- field theory. In the next section [3] we consider the BCS-type mean- 
field theory in homogeneous space and describe the effect of boundary conditions on the 
stability of different phases. In section H] we consider the effects of trapping potential 
using the BdG method. In section [5] we show the main results obtained from the 
BdG method, especially regarding the coexistence of two pairing gaps. We conclude by 
discussion in section |6j 
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2. The system setup 

The general mean-field Hamiltonian for a three-component system in the contact 
interaction potential approximation is (up to a constant) 

h 2 v 2 



H MF = I d^liv) 

~ — 1 O Q J 



ct=1,2,3 



2m a 



V a {v)-^ + W a {v) 



+ l -Y, I d 3 r A aa/ (r)^i(r)¥,(r) + h.c, (2) 

where the first term of the Hamiltonian includes contributions from the kinetic energy, 
the external trapping potential V a {v) (that can depend on the component |cr)), and 
the chemical potentials /i CT , respectively Interactions are described by the two-body 
scattering T-matrix. In the contact interaction potential approximation it can be written 

as 

^ry) = ^^(r-r'), (3) 

where a aa > is the scattering length between atoms in hyperfine states \a) and \a'), and 
m r = 1ra a m a i / {m a + m a r) is twice the reduced mass. The Hartree fields are denoted by 
W a {Y) = J2a^a> Kr(7'( r ) n (7'( r ) an d the densities are n a {r) = ($f^(r)ty a (r)) . The pairing 
(mean-)field A CTCT /(r) = f/ crcr '(r)(\l/ -(r)^' (T /(r)) includes a renormalized interaction U aa >(r) 
that is used to remove the ultraviolet divergence following the standard procedure (see 
below). In our model we neglect the possibility of three-body bound states and other 
three-body effects that can affect the lifetime of the gas [7j. 

A three-component system has three possible pairing fields corresponding to the 
three interaction channels Uu, Ui$, U23, and these can be combined into a total pairing 
field vector A = [A23, — A 13 , A 12 ] T . Identical properties make the system SU(3)- 
symmetric and the pairing vector can be reduced to a single gap by a simple unitary 
transformation. The orientation of the pairing vector corresponds to a choice of the 
global gauge [HI [T7J [HI [191 125] , and the simplest choice is the one where only one of 
the pairing gaps, say A 12, is nonzero. This of course makes atoms in component |3) 
effectively noninteracting. Indeed, in an SU(3)-symmetric case, there always exists a 
gapless branch describing unpaired atoms. 

In our study one interaction is always suppressed (we choose C/13 = 0). This 
is the case for example in 6 Li where Feshbach resonances between the three lowest 
hyperfine states |1) - |2) (at B = 834 G) and \2) - |3) (at B = 811 G) lie close to 
each other, while the resonance |1) — |3) (nearest one at B — 690 G) is sufficiently 
far away [27]. Similar behavior occurs in 40 K [28], where the richness of the hyperfine 
level structure allows even more freedom in choosing the suitable interaction strengths. 
Moreover, mixtures of 6 Li and 40 K offer interesting possibilities [29[[30]. Also we do not 
consider here interactions between |1) and |3) induced by the component |2) [31]. Thus, 
neglecting the |1) — 13) interaction channel altogether, the symmetry is broken at least to 
SU(2) xSU(l). Analogously to the SU(3) symmetric case, the total pairing field is now a 
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two-dimensional vector A = [A 12 , A 23 ] T and in the symmetric case, where components 
|1) and |3) are identical, it is preserved under spin rotations of the hyperfme states |1) 
and 1 3). This symmetry implies that the ground state is degenerate with respect to 
the orientation of the total pairing field vector. However, the degeneracy is lifted by 
changing masses, chemical potentials or interaction strengths and, as we will soon show, 
also by imposing boundary conditions such as fixing the number of atoms in different 
components. 

Boundary conditions, such as fixed particle numbers or fixed chemical potentials, 
manifest themselves in different ways in atomic gases. While the total particle numbers 
are, in practice, fixed, the local densities are not as the particles are allowed to move 
around in the trap. Hence, from the local density approximation (LDA) point of 
view, locally the relevant boundary condition appears to be a fixed chemical potential. 
However, globally the relevant boundary condition is the fixed particle number, and in 
the BdG method we indeed fix the mean particle number. Below we will also discuss 
how the two pictures merge in the limit of large system size iV — > oo. 



3. Homogeneous system 



For a homogeneous system the densities and gaps lose their spatial dependence. This 
corresponds to the usual BCS-approximation in which pairing can only occur between 
atoms with opposite momenta |k, a) and | — k, a'). The mean- field Hamiltonian can be 



written in matrix form as 

T 
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(4) 



3k 



/ 



where C is constant and the single-particle dispersion is £ CTk = |^ — We have 
here neglected the Hartree fields since at the level of our approximation they provide 
only a constant energy shift in a homogeneous system. The standard approach calls for 
diagonalizing this using the Bogoliubov transformation, and iteratively solving for the 
pairing fields A 12 and A 23 . In order to satisfy the fixed mean atom number boundary 
condition, the iteration must adjust the chemical potentials in a self-consistent manner 
as well. Notice that the inherent atom number fluctuations implied by the mean-field 
theory play no role here as long as the typical fluctuations (scaling as y/~N) are much 
smaller than the atom numbers in different species (scaling as N). 

As discussed above, in the symmetric case = /x 3 , m 1 = m 3 , U 12 = U23), the 
Hamiltonian has SU(2)xSU(l) symmetry and all the pairing fields with A^ 2 + A| 3 
constant yield the same total energy. The ground state is thus degenerate. However, 
different orientations of the pairing field vector yield different atom numbers in 
components |1) and |3). Thus, fixing the numbers of atoms N% and iV 3 breaks the 
degeneracy and a well-defined energy minimum is found. Figure [1] shows typical energy 
landscapes (Hmf + J2a /V^o-) as a function of the pairing fields Ai 2 and A 23 for equal 
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A 12 [E F ] A 12 [E F ] 

Figure 1. Logarithmic energy landscape with constant particle numbers (a) Ni — 

N 3 = 1.0iV 2 (b) Ni = N 3 = l.2N 2 (c) N x = N 2 ,N 3 = 1.1N 2 (d) Ni = 1AN 2 ,N 3 = 

1.2N 2 . The interaction strengths are equal (/cFa^) -1 = (kpa^)^ 1 = —0.5, (kp = 
(67r 2 n 2 ) 1/3 ). 



interaction strengths Uu = U^- The Fermi momentum kp here and throughout this 
work is defined as the Fermi momentum of the component |2), kp = y^2mEp =2 / h, where 
Ep is the Fermi energy of the component \o~). The energies have been calculated for 
fixed atom numbers: the chemical potentials \i a are solved for every point (Ai 2 , A23) so 
that the atom number constraints are satisfied. The figures show clearly how the ground 
state becomes non-degenerate and realizes itself in a particular combination of pairing 
fields. In the special case where there is equal number of atoms in all three components, 
Ni = N 2 = N 3 , the ground state is still degenerate. However, if the number of atoms in 
component |2) is changed (keeping N\ = N% but N2 7^ N±), the degeneracy is broken and 
a non-degenerate energy minimum appears. This is in stark contrast to the case where 
the chemical potentials are kept constant and the atom numbers are allowed to vary. In 
such a case the ground state remains degenerate as long as the chemical potentials for 
components |1) and |3) are equal, Hi = fi 3 . 

In the case of a number mismatch or difference in the interaction strengths of the 
components |1) and |3), the energy minimum will be shifted from the equal pairing case. 
Depending on the number of atoms in component |2), the minimum appears either at the 
edge of the energy landscape (yielding either of the two pairing fields A12 or A23 zero) or 
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Figure 2. The single-particle occupation numbers for Nx = N3 = I.2N2 and equal 
pairing fields A12 = A23 = 0.25-Ef- All atoms in the hyperfme state |2) are paired but 
part of the atoms in components |1) and |3) are unpaired. The unpaired atoms form 
a well-defined Fermi sphere resulting in a step at the Fermi surface. 



somewhere in between. This too is an important difference to the case of fixed chemical 
potentials where the breaking of the symmetry (by either changing chemical potentials 
or interaction strengths) always results in either of the two pairing fields dominating 
and the other becoming zero. Thus, for fixed chemical potentials one does not observe 
coexistence of the two pairing fields A i2 and A 23 except possibly in the symmetric, or 
degenerate, case, whereas for fixed atom numbers the coexistence phase (described by 
two non-vanishing pairing field order parameters A 12 and A23) is very real. 

The pairing scheme is revealed by the momentum distribution of each state. 
Figures |2] and [3] show the momentum distributions of the three components for equal 
pairing gaps (Ai 2 = A 23 ) and for projected pairing gaps (obtained by setting A 23 = 0, 
which is always an allowed solution, and minimizing the energy by varying only A 12 .) In 
the first case, equal pairing gaps imply that there are equal numbers of 12 and 23 pairs. 
Since only zero-momentum Cooper pairs are considered here, one can filter the paired 
atoms from the momentum distributions and determine the momentum distribution 
of unpaired atoms by calculating the difference (n a k) — (^2fc)/2 for a = 1,3. The 
distribution of unpaired atoms is seen to form a clear Fermi sphere, but with maximal 
occupation probability of 0.5. In the case of projected pairing gap in Figure |3] the 
pairing atoms |1) and |2) form a breach due to a number mismatch between the two 
components. 

The boundary condition of fixed atom numbers for each component separately is 
natural for atomic gas experiments. However, the results for a homogeneous gas must 
be approached with caution since experiments are always conducted in nonuniform 
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Figure 3. Single-particle occupation numbers as in Figure [2] but now for projected 
pairing field vector A23 = 0, A12 = 0.32-Ep. All atoms in the component |3) are 
unpaired, revealing a noninteracting Fermi sea. Due to the number mismatch between 
components |1) and |2), excess atoms in |1) will form a breach. 



trapping potentials. Using these homogeneous system results in conjunction with local 
density approximation means locally fixing the chemical potentials instead of the atom 
numbers. This discrepancy on which boundary condition to use can be solved by treating 
the trapping effects explicitly using the Bogoliubov-deGennes method. 

4. Harmonic trap — the Bogoliubov-deGennes method 

In order to consider trapped systems, we use the Bogoliubov-deGennes method that 
allows the inclusion of trap effects exactly. The mean-field BdG method is not expected 
to be able to capture all relevant physics in the strongly interacting regime. However, 
in an imbalanced two-component system, it has been shown [52] that, for small 
polarizations and symmetric trap geometries, there is a good agreement between the 
mean-field BdG approach and real-space dynamical mean-field theory. We solve the 
three-component mean- field system in a spherically harmonic trap V a (r) = \m a bj 2 a r 2 
using the eigenbasis of the 3-dimensional harmonic oscillator 

^ a (r) = J2Kn(r)Y lm (n)c nlma , (5) 

nlm 

where Y\ m are the spherical harmonics and the radial wavefunctions are given by 

Kn(r) = ^ m ^T ,A \[^ ( 6 ) 
Here Li +1 ^ 2 (f^) is the associated Laguerre polynomial and f a = r^/m^uj^/h. 
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The mean-field Hamiltonian separates for different /-quantum numbers H = 
^[ Hi + C, such that [H, Hi] = and C is constant. Introducing a finite cutoff energy 
E c and keeping only single-particle states with energy less than the cutoff allows writing 
each Hi using block matrices 



H 



( c f \ 

\4ij 



I 6/1+ J 



12 



12 



12 



-12 



23 







23 

E /3 + ^3; 



\ 
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c 2« 
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V C 3i 
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(7) 



where Jo 



J l 21 + JJ,,. The block matrices are defined as 
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/ J! 
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NO 



J 1 



Jl 
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'NN 
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J 



( F l 



and Fl 



NO 



'ON 



'NN 



J 



with the Hartree shift 



poo 

J l M = U*, / drr 2 B? an (r)n a ,(r)KAr) 
Jo 

and the pairing field 



alrr 2 R l an {r)K„'{r)R l a , n ,{r) 



(8) 



(9) 



(10) 



The connection between the interaction strength U^ a , used in the Hartree shift and 
the bare interaction strength U aa > will be discussed below. The energy matrix e la is 
diagonal with elements e anl = ftw a (2n + I + 3/2) — \i a and the operator vectors are 
c al = [c a0l0 • ■ ■ c . ivro ] T . We denote the number of single-particle states with fixed I, 
whose energy is below the cutoff, as Ni = [E c j (hu) —1 — 3/2} /2. 

Similarly to free space, the Hi matrices can be diagonalized using the Bogoliubov 
transformation which is provided by unitary 3N l c x 3^-matrices W ; . By inserting 
the identity operator (W')tW z between the matrix and the operator vectors, we have 

. The rotation matrix W ; is 



W z 



C ll C 2l C 3l 



r 



the quasiparticle basis as ( j u y 2l 7 3Z 
chosen such that the matrix in is diagonalized. 

The equations for the pairing fields and the densities are 

A CTCT ,(r) = U aa , ^(2/ + l)RUr)K' n '(r)J2 W U + n,M'N l+ n'M E 
and 



Ji- 



nn' l 



n a (r) = £(2i + l)R l an (r)R l an ,(r) £ Wi Ni+nj Wi Ni+n , ^n F {{-lf E 3 ) 



11) 



(12) 



where n-p is the Fermi distribution and a — a — 1, a £ {1, 2, 3}. The total number of 
particles in each component \a) is obtained by integration N a = J °° dr r 2 n a {r). 

As in usual BCS theory, the gap equation is ultraviolet divergent, hence the energy 
cutoff E c . In order to make the model cutoff independent, we follow a standard 
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approach [33] and use a renormalized interaction U aa i but now generalized to a three- 
component system 

1 1 m r k c , aa >(r) 



-cw(r) (13) 



where 



Here the momentum cutoff k CjCra i and local Fermi momentum k-p^a' are defined as 

h 2 k 2 r ,(r) h 2 B,(r) 

2 Zr [Ec ~ KAr)] ^ 2m r = i 2 ^Ar) ~ W a (r) - WAr)] , (15) 

where the local average chemical potential is 

fl aaf (r) = W - V„(t) + n a , - VAr)] /2. (16) 

We solve these gap and number equations self-consistently using fixed point 
iteration. For every iteration step in the gap equation, we solve the chemical potentials 
fi a to keep the particle numbers constants. The iteration is terminated when the 
subsequent gap profiles in the iteration differ by at most 5 x 10~ 5 fad. The cutoff energy 
is chosen to be 2.5 x ma.x{Ep} , (with a higher cutoff, the results do not qualitatively 
change). If the convergence is slow, we try different initial values to ensure that the 
final result is correct. We use a small finite temperature (T = 10~ 3 Tp) to smoothen the 
Fermi distribution and to help solving the number equations in presence of a discrete 
energy spectrum. However, we have checked that the results are unchanged even for zero 
temperature for example in Figures [6] and [7J We do not consider higher temperatures in 
this work, however, we have checked that our results are sufficiently robust to survive 
low but experimentally relevant temperature T = 0.05 Tp. An example of the effect 
of the temperature is shown in Figure H] where exactly the same parameters as in 
Figure [6] are used except that T = 0.05 Tp. The results, especially the coexistence 
region, are practically identical, only the minor features at the edge of the gas have 
been smoothened. 

The Hartree fields become infinite with a diverging scattering length a acr / — > oo. 
This unphysical effect is caused by improper treatment of two-body scattering effects 
and in practice these energy shifts are limited by the Fermi energy. Monte Carlo results 
on a two-component Fermi gas suggest that the Hartree fields at unitarity do not exceed 
\W\ ~ 0.5^ [Ml EH]. We limit the Hartree field interaction to be smaller than this 
value by imposing a hard cutoff on the Hartree interaction strength. That is, instead of 
the bare interaction XJ aa i we use U^., which is limited from above by 

I^MO)! <0.5<. (17) 

Notice that the component |2) experiences two Hartree fields due to the two components 
|1) and |3) and thus the total Hartree shift experienced by this component can be 
up to double the above cutoff. Since the Hartree shifts induce a mismatch between 
the Fermi surfaces, the pairing amplitude is reduced. Notice that this is in contrast 
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Radial distance r[a usc ] 

Figure 4. Typical gap and density profiles in harmonic trap calculated for temperature 
T = 0.05 Tp. Other parameters are the same as in Figure [6] 

to the balanced two-component case in which both components experience the same 
Hartree potential and the densities remain thus equal. In three component systems, 
the inclusion of the third interaction {7 13 would reduce the mismatch, but because it is 
usually weaker, the mismatch does not totally disappear. However, the mismatch can 
be countered by careful choice of interaction strengths and atom numbers, so that local 
density imbalances are reduced. 

The Hartree shift can be seen to produce interesting shell structures for certain 
parameters. Such exotic shell structures created by the Hartree shift will be considered 
elsewhere, requiring a more complete treatment of the Hartree effect to confirm their 
validity. In the present work, we have chosen to focus on a parameter range in which 
these peculiar features are not present and in all of the results shown in this work except 
in Figure [5] we neglect the Hartree effects, i.e. we use U^ a , = for all a, a'. Our choice 
does not significantly limit the parameter range, because these effects appeared always 
only in tiny islands in the parameter space, typically at the edges of the trap. Indeed, 
we have checked that the inclusion of the Hartree shift in the way described above does 
not qualitatively change the results presented here. An example is shown in Figure El 
which presents the case of Figure [6] but with Hartree fields included. The qualitative 
behaviour is the same although the numerical values of the order parameters are smaller. 
Importantly, the Hartree fields do not affect the coexistence of the two order parameters. 
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Figure 5. Typical gap and density profiles in harmonic trap calculated for zero 
temperature but including the Hartree energy shift in the way described in the main 
text. The parameters are the same as in Figure [6] 




Radial distance r[a osc ] 



Figure 6. Typical gap and density profiles in harmonic trap showing a large 
coexistence region of the pairing fields A12 and A23. The parameters are N 2 — 3 x 10 4 , 
Ni = 0.8N 2 , N 3 = 0.7N 2 , and {k F a^)" 1 = (fc F a 23 ) _1 = -0.50. 
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Radial distance r[a osc ] 

Figure 7. Gap and density profiles as in Figure H] but for mismatched interaction 
strengths (^012) 1 — —0.50, (kFa2z)~ 1 = —0.45. The two pairing fields A12 and A23 
separate into different regions, with a A23 core surrounded by a A12 shell. 

5. Results 

Figures [6] and [7J show typical gap and density profiles obtained from the BdG method. 
In the former the gaps follow the density distributions, and the two pairing gaps Ai 2 
and A 23 are present across the trap. In the latter the gaps are spatially separated, with 
the A12 pairing field concentrated at the edge of the trap and A 2 3 at the center of the 
trap. There is no clear interface between the two pairing regions, but the penetration 
length of pairing field A 12 inside A 23 is relatively constant when increasing the system 
size (increasing atom numbers N), characteristic length scale given by the oscillator 
length r osc = ^Jh/ni2U)2- Locally the dominating pairing channel is the one for which 
the atom densities are least mismatched, the strength of the interaction being only a 
secondary factor. This local nature of pairing allows interesting shell structures [36] as 
shown in Figure [7J However, we do not pursue these issues here but rather concentrate 
on more general features. Notice that the number of atoms in component |3) has been 
chosen to be slightly smaller than in components |1) and |2), but the features shown 
here are very general. 

Figures [HI a) and b) show the pairing fields Ai 2 and A 23 as a function of interaction 
strength (/cpo^) -1 . When either of the two interaction strengths ai 2 , a 23 is significantly 
stronger, the corresponding pairing channel will dominate. The crossover between the 
two regions occurs at (£;Fa 23 ) _1 = —0.46 (not at (/cfg^) -1 = —0.50 because of the atom 
number mismatch Ni > iV 3 ). To better characterize the coexistence of the two pairing 
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Figure 8. (a) Pairing field A12 (in units of hut) as a function of the interaction strength 
U23- (b) Pairing field A23. (c) The coexistence parameter P co — A 12 A 23 /(Ai 2 + A 23 ) 
shows the coexistence areas. Atom numbers in a), b) and c) are the same as in Figure|6] 
In (d) the coexistence parameter is plotted at r = for different numbers of particles 
N. 



gaps, we define a dimensionless coexistence parameter 

A12A23 



A 2 



A 2 ' 

^23 



(18) 



Figure |8] c) shows this parameter as a function of interaction strength and position, 
revealing a large coexistence region in the somewhat narrow interaction strength window 
—0.51 < (k F a 2 3)~ 1 < —0.47, but also a coexistence region close to the edge of the 
trap across a wide range of interactions. In Figure |8] d) we show how the coexistence 
parameter at the center of the trap scales with increasing system size N (the atom 
numbers Ni and iV 3 are scaled correspondingly so that the relative polarizations are 
fixed). The coexistence area is suppressed as N grows large, implying that the 
coexistence may vanish in the thermodynamic iV — > 00 limit. However, with sufficiently 
accurate choice of interaction strengths, the coexistence region should be experimentally 
accessible with reasonably sized atom gases. We have not studied the scaling of the 
coexistence regions at the edge of the trap, but since the penetration length in Figured 
is given by the oscillator length we expect the N — > 00 limit to yield a phase separation 
into a A 2 3 core and a surrounding A 12 shell. 

To better understand the nature of the pairing scheme in the coexistence areas, 
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Figure 9. Occupation numbers (n CTJi ,; = o) = (c an0 c an o) for each component a as a 
function of the quantum number n with I = 0. The parameters have been chosen 
symmetrically Ni = N 3 = 16000, N 2 = 20000, and (fea^)" 1 = (k F a 23 y 1 = -0.50, 
resulting in identical pairing gaps Ai 2 (r) = A 2 3(r). 

Figure [9] a) shows the occupation numbers of different n-quantum number states for 
a symmetric case A 12 (r) = A 13 (r). The angular momentum quantum number / = 
chosen here acts as a representative of a more general behavior for general I. The figure 
reveals a very similar pairing scheme as in the homogeneous system, see Figure [2j The 
unpaired atoms are distributed among the two components |1) and |3) and form a step 
at the Fermi surface. 

6. Conclusions 

We have studied the pairing of a three-component Fermi gas when two of the interspecies 
interaction channels are dominant. In a homogeneous system, we showed that different 
boundary conditions, namely fixing the chemical potential or fixing the atom numbers, 
produce qualitatively different results and phases. We have not considered the possibility 
of a phase separation in which the two pairing fields would be spatially separated in 
otherwise uniform system. 

For trapped systems, our BdG study reveals an interesting coexistence region where 
both pairing channels A 12 and A 23 are present. This is a mesoscopic effect and likely 
to vanish in the limit of a large system N — > oo resulting in phase separation into shells 
of different pairing fields. However, the coexistence region is present at atom numbers 
relevant for atom gas experiments, making the observation of this intriguing double-gap 
prediction feasible. 

There is already a wide range of standard experimental techniques to detect such 
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pairing correlations. For example, one could use radio-frequency spectroscopy [371 EH 
[39] for driving atoms from the hyperfine states |1) and |3) separately into some fourth 
noninteracting state |e). Other possibilities include transforming pairing correlations 
into molecular 12 and/or 23 pairs through magnetic field sweeps [40] or optical molecular 
spectroscopy jH]. 
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